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IT} ; ABSTRACT 

We will discuss here an indirect lattice evaluation of the baryon axial singlet current 
matrix element. This quantity may be related to the fraction of nucleon spin carried 
by the quarks. The appropriate structure function has recently been measured 
(EMC experiment). As in this experiment, we hnd that the quarks do not appear 
to carry a large component of the nucleon spin. 
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1. Introduction and Theoretical Discussion 



Hadrons appear to be far more complicated than the (rather successful) con- 
stituent quark model would suggest. For example an old result is from the nN 
sigma term, which seems to give a rather large strange component to the nucleon 
mass. A more recent result is from the EMC experiments which suggests that 
constituent quarks are responsible for very little of the nucleon spin. These are 
non-perturbative effects and so is an area where lattice calculations may be of some 
help. 

The EMC experiment measured deep inelastic scattering (DIS) using a /i beam 
on a proton target. The new element in the experiment was that both the inital \i 
and proton were longitudinally polarised. A measurement of the difference in the 
cross sections for parallel and anti-parallel polarised protons enabled the structure 
function g\(x. s Q 2 ) to be found. Theoretically! this is of interest as using the Wilson 
operator product expansion, moments of the structure function are related to certain 
matrix elements. In this case, defining s^Aq = (Ps| 57^75(7 |Ps) for q = u, d or s 
quarks (i.e. the expectation value of the axial singlet current) the lowest moment 
is then given by 
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« 0.20 + -As, 



(1) 



(Q 2 emc ~ HGeV 2 ). Results from neutron and hyperon decays have been used to 
eliminate two of the unknowns on the RHS of this equation. 

Aq can be given a physical interpretation, as the quark spin operator, Sf uark , 
is local and gauge invariant [in distinction to the gluon spin operator or orbital 
angular momentum operator)tl and leads to 

2(Ps\Sl uark \Ps) = AS = Au + Ad + As 

« 0.68 + 3As. (2) 

Before the EMC experiment, the Ellis-Jaffe sum rule was applied: the strange 
content of the proton was taken as zero, i.e. As = 0. However the EMC experiment 
gave for the LHS of eq. ([[]) ~ 0.126 giving As « —0.2, a large negative strange 
contribution to the proton and S « 0, i.e. the total quark spin is zero. This 
situation is often referred to as "the spin crisis of the quark model" , although what 
was actually measured was the nucleon matrix element of the axial singlet current. 

2. The lattice calculation 

We now turn to our lattice calculation. By working in euclidean space on a 
lattice we can turn our problem into a statistical mechanics one, which can be 
approached numerically via the calculation of correlation functions using Monte 
Carlo techniques. We have generated configurations^ using dynamical staggered 
fermions (x) on a 16 3 x 24 lattice at j3 — 5.35, m = 0.01, which corresponds to 
a quark mass of about 35MeV. Staggered fermions describe 4 degenerate quark 
flavours, so we do not have quite the physical situation that we wish to describe 
(i.e. 2 light quarks and one heavier quark). Nevertheless this discretisation of the 
fermions has some advantages - most notably good chiral properties (as m — > 0) 
when the chiral symmetry is spontaneously broken and there is a Goldstone boson, 
the 7r. 

The general procedure on the lattice to measure matrix elements is to consider 
3-point correlation functions!. For an arbitrary operator Q, it can then be shown 
that (B is the nucleon operator) 

C(t;r) = <£fe*(t)n(r)fl^,(0)> 

= J2 A aP (a\m^ T ^, (3) 

a,(3=N,K 

for \T ^> t ^> t 3> 0, where T is the time box size = 24 here, /i^ = +exp (—E N ) 
and the parity partner (which always occurs when using staggered fermions) is 
denoted by A, so that fi\ = — exp (— E\). The amplitudes A a p and energies E a 
are known from 2-point correlation functions. Thus from eq. (|3|) we can extract 
(Ps\Cl\Ps). Simply setting fl to be X7i7sX nas 5 however, a number of disadvantages: 



the correlation function in eq. (|3|) has connected and disconnected parts, in the fit 
there are cross terms present, and the operator must be renormalised. Another 
approach is to consider the divergence of the current^. This can be related to the 
QCD anomaly. In the chiral limit an equivalent formulation of the problem is thus 

(n f = 4) 

AS = Urn \ ^ P(p),s\Q(p-p')\P(p'),s'}, (4) 
p'^p i{p — p') ■ s 

with topological charge 

Q(p) = / d 3 xF^e^ s . (5) 

(We shall take p' = 0, p = (0,0,p).) On the lattice we have used the Liischeri 
construction for this charge. Although technically complicated it is integer valued 
and hence has no renormalisation. In the 3-point correlation function there are no 
cross terms (Q has a definite parity). A disadvantage is that we need to take p — > 0. 
On our lattice the smallest momentum available is the rather large p ~ 500MeV. 

Nevertheless on attempting this measurement, we find a reasonable signal with 
AE rs 0.18(2). This is small and tends to support the EMC result, which would 
indicate a rather large sea contribution to the proton. The existence of the QCD 
anomaly also proved important. Further details of our calculation are given inS. 
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